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ABSTRACT 

We interpret recent ALMA observations of 2 > 6 normal star forming galaxies by 
means of a semi-numerical method, which couples the output of a cosmological hydro- 
dynamical simulation with a chemical evolution model which accounts for the con¬ 
tribution to dust enrichment from supernovae, asymptotic giant branch stars and 
grain growth in the interstellar medium. We find that while stellar sources dominate 
the dust mass of small galaxies, the higher level of metal enrichment experienced by 
galaxies with M sta r > 10 9 M 0 allows efficient grain growth, which provides the domi¬ 
nant contribution to the dust mass. Even assuming maximally efficient supernova dust 
production, the observed dust mass of the z = 7.5 galaxy A1689-zDl requires very 
efficient grain growth. This, in turn, implies that in this galaxy the average density 
of the cold and dense gas, where grain growth occurs, is comparable to that inferred 
from observations of QSO host galaxies at similar redshifts. Although plausible, the 
upper limits on the dust continuum emission of galaxies at 6.5 < z < 7.5 show that 
these conditions must not apply to the bulk of the high redsliift galaxy population. 

Key words: dust, extinction galaxies: evolution galaxies: higli-redshift galaxies: 
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1 INTRODUCTION 

Observations at millimeter (mm) and sub-millimeter (sub- 
mm) wavelengths have provided convincing evidence of 
rapid dust enrichment at high redshift (for a recent re¬ 
view see Carilli & Walter 2013 and references therein). Dust 
masses as large as ~ 10 s M@ have been inferred for 5 < z < 7 
QSO host galaxies, requiring a very efficient dust formation 
channel that must operate in less than 1 Gyr of cosmic evo¬ 
lution. Theoretical studies have shown that although stel¬ 
lar sources of dust, including both Supernovae (SN) and 
Asymptotic Giant Branch (AGB) stars, can be fast enough 
to significantly contribute to high redshift dust enrichment 
(Valiante et al. 2009), the observed dust masses require ef¬ 
ficient grain growth in the dense phase of the interstellar 
medium (ISM, Michalowski et al. 2010; Valiante et al. 2011; 
2014). Until very recently, sub-mm/mm continuum observa¬ 
tions at 2 ~ 6 — 7 have been limited to “extreme” galax¬ 
ies, QSO hosts or strong starbursts, which are characterized 
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by star formation rates ranging from several hundreds to 
thousands of solar masses per year. However, at these high 
redshifts we expect the bulk of the galaxy population to 
be characterized by “normal” galaxies, which form stars at 
rates from a few to a few tens of solar masses per year. Deep 
ALMA and Plateau de Bure Interferometer (PdBI) follow 
up observations have mostly provided upper limits on the 
dust mass of such galaxies (Kanekar et al. 2013; Ouchi et 
al. 2013; Ota et al. 2014; Schaerer et al. 2015; Maiolino et 
al. 2015), with the exception of the gravitationally lensed 
a = 7.5 Lyman Break Galaxy (LBG) A1689-zDl, whose 
dust continuum emission has been recently detected with 
ALMA (Watson et al. 2015). The estimated dust mass is 
4±| x 10' Mq , showing that dusty galaxies have already 
formed at z > 7. 

I11 this paper, we investigate early dust enrichment in 
“normal” star forming galaxies at 2 > 6. We use the out¬ 
put of a cosmological hydro-dynamical simulation which al¬ 
lows the prediction of the gas, stellar and metal content of 
galaxies along their hierarchical assembly (Maio et al. 2010). 
We then post-process the simulation output with a semi- 
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Table 1. Physical properties of the galaxy sample collected from 
the literature. Lensed objects have been corrected using the ap¬ 
propriate magnification factor (// = 4.5, 9 and 9.3 for HCM6A, 
A1703-zDl and A1689-zDl, respectively). All the quantities are 
computed assuming conversion factors, scaling relations and dust 
properties presented in the text. References: a: Schaerer et al. 
(2015); b: Kanekar et al. (2013); c: Ota et al. (2014); d: Ouchi et 
al. (2013); e: Maiolino et al. (2015); /: Watson et al. (2015). 


Name 

z 

M uv 

Log Mstar 

Log M dus t 



mag 

[Mq] 

[Mq] 

A1703-zDl a 

6.800 

- 20.3 

9.2 ±0.3 

< 7.36 

z8-GND-5296“ 

7.508 

-21.4 

9.7 ±0.3 

< 8.28 

HCM6A 6 

6.560 

-20.8 

9.5 ±0.3 

< 7.61 

IOK-l c 

6.960 

-21.3 

9.7 ±0.3 

< 7.43 

Himiko d 

6.595 

-21.7 

9.9 ±0.3 

< 7.30 

BDF-3299 e 

7.109 

-20.44 

9.3 ±0.3 

< 7.02 

BDF-512 e 

7.008 

-20.49 

9.3 ±0.3 

< 7.36 

SDF-46975 e 

6.844 

-21.49 

9.8 ±0.3 

< 7.38 

A1689-zDll 

7.500 

-19.7 

9.0 ±0.3 

7.51 ±0.2 


analytical chemical evolution model with dust (Valiante et 
al. 2009; 2011; 2014; de Bennassuti et al. 2014) to estimate 
their dust masses. A similar semi-numerical approach has 
been followed by Dayal, Hirashita & Ferrara (2010), who 
considered the contribution to dust enrichment of SN and es¬ 
timate the Far InfraRed (FIR) detectability of liigh-redshift 
galaxies classified as Lyman-a emitters (LAE). More re¬ 
cently, Hirashita et al. (2014) have been able to put con¬ 
straints on the dust production rate by SN using the upper 
limit on the dust continuum emission of Himiko, one of the 
best-studied LAE at 2 ~ 6.6 (Ouchi et al. 2013). Here we 
attempt to improve on these previous studies by considering 
dust enrichment from SN, AGB stars and grain growth in 
the ISM with the aim of assessing their relative contribution 
to the total mass of dust in the ISM of z > 6 galaxies. We 
compare the model predictions with currently available up¬ 
per limits on the dust mass at 6 < z < 7.5 and discuss the 
implications posed on the models by the newly discovered 
dusty galaxy at z = 7.5 (Watson et al. 2015). 

The paper is organized as follows: in section [2] we intro¬ 
duce the sample of observed galaxies that we have collected 
from literature papers; in section [3] we give a brief presenta¬ 
tion of the semi-numerical model and in section[4]we discuss 
the main results. Finally, in section[5]we draw our main con¬ 
clusions. 


2 THE OBSERVED SAMPLE 

In Tabic |T] we have summarized the properties of the sam¬ 
ple collected from the literature. For 8 out of the 9 galax¬ 
ies, upper limits on the dust continuum emission have been 
obtained by means of deep ALMA and PdBI observations 
(Kanekaer et al. 2013; Ouchi et al. 2013; Ota et al. 2014; 
Schaerer et al. 2015; Maiolino et al. 2015). The z = 7.5 
LBG A1689-zDl is currently the most distant UV-selected 
galaxy for which a dust continuum detection has been ob¬ 
tained with ALMA (Watson et al. 2015). For all these galax¬ 
ies, the star formation rates estimated from the UV luminos- 


iuQ range between 4.1 and 22.4M©/yr. Following Schaerer 
et al. (2015), we derive the stellar masses using the mean 
relation between the UV magnitude and the stellar mass ob¬ 
tained by the same authors from detailed fits of the spectral 
energy distribution of a sample of z ~ 6.7 LBGs, including 
nebular emission and dust attenuation: 

Log(M st ar/M©) = -0.45 x (Muv + 20) + 9.11. (1) 

Dayal et al. (2014) obtain a similar M sta r — Muv relation on 
the basis of a theoretical model, which also accounts for the 
rapid decline of the UV luminosity with stellar ages. 

Assuming that the dust is optically thin in the rest- 
frame FIR, the upper limits on the dust mass have been 
obtained taking into account the effect of the Cosmic Mi¬ 
crowave Background (CMB) on the intrinsic dust emission 
(da Cunha et al. 2013). In Table [T] we report dust masses 
adopting a dust temperature^ Td us t,o = 35 K and a dust 
emissivity k„ rea = ko(Ao/A res ) 0 with ko = 0.77cm 2 /gr, 
Ao = 850 fi m and /3 = 1.5 (Ota et al. 2014). The dust mass 
can increase/decrease by ~ 0.4 dex for variations in the dust 
temperature in the range 25K < Td ua t,o < 45K. In addition, 
the poorly constrained dust properties make the value of 
the emissivity coefficient to be adopted very uncertain (see 
Table 1 in Hirashita et al. 2014). At fixed dust tempera¬ 
ture, we find that variations in k^ among values adopted 
in the literature (Weingartner & Draine 2001) or applied to 
submm observations of high-z galaxies (Michalowski et al. 
2010; Valiante et al. 2014; Watson et al. 2015) introduce an 
additional ~ 0.3 dex uncertainty in the estimated dust mass. 


3 THE MODEL 

We briefly describe the simulation used in this work and 
interested readers are referred to Maio et al. (2010), Camp- 
isi et al. (2011), Salvaterra et al. (2013) and Dayal et al. 
(2013) for more details. The smoothed particle hydrody¬ 
namic (SPH) simulation is carried out using the TreePM- 
SPH code GADGET-2 (Springel 2003) and has a periodic co¬ 
moving box size of 30 h~ 1 cMpc. It initially contains 2 x 320 3 
particles, with the masses of gas and DM particles being 
9 x 10 6 1i _1 Mq and 6 x 10Ti“ 1 Mq, respective!^. The code 
includes the cosmological evolution of both atomic and ion¬ 
ized Hydrogen, Helium and Deuterium (e.g. Yoshida et al. 
2003; Maio et al. 2009), PopIII and PopII/I star formation 
according to the corresponding initial mass function (IMF; 
Tornatore et al. 2007), gas cooling from resonant and fine- 
structure lines (Maio et al. 2007) and feedback effects as¬ 
suming a multi-phase ISM (Springel & Hernquist 2003). The 

1 We have adopted a conversion factor of /Cuv = 1.15 X 
10 -28 Mq yr —1 /ergs —1 Hz _1 , as appropriate for a stellar popula¬ 
tion with metallicity 0.3 Zq < Z < Zq formed at a constant SFR 
with a Salpeter IMF in the range [0.1 —100] Mq at age ~ 300 Myr 
(Madau & Dickinson 2014). 

2 Tdust 0 I s the dust temperature heated by a stellar radiation 
field at z = 0 (da Cunha et al. 2013). 

3 The cosmological model corresponds to the ACDM Universe 
with DM, dark energy and baryonic density parameter values of 
(f2 m , 12 a, Dt,) = (0.3,0.7,0.04), a Hubble constant Ho = 100/r = 
70km s — tMpc^ 1 , a primordial spectral index n s = l and a spec¬ 
tral normalisation ers = 0.9. 
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Figure 1. Predicted dust masses of the simulated galaxies as 
a function of the stellar mass. For each galaxy, the dust mass 
without (with) grain growth is shown by a square grey (circle 
blue) point (see text). The adopted grain growth timescale is 
Tacc.o = 2Myr. In the lower panels, the reverse shock destruction 
of SN dust is neglected. For the sake of comparison, we have 
reported the same data points shown in Table 1 in the two panels: 
Schaerer et al. (2015, squares), Maiolino et al. (2015, triangles) 
and Watson et al. (2015, circle point). 


Figure 2. Redshift evolution of the dust mass for the simulated 
galaxies with stellar masses in the range Log M s t a r/Mg > 9. Each 
line represents the average contribution of all the galaxies with the 
shaded area indicating the dispersion among different evolution¬ 
ary histories. The lower, intermediate and upper lines show the 
contribution to the total mass of dust of AGB stars, stellar sources 
and grain growth with an accretion timescale of r aC c,o = 2 Myr 
(upper panel) and 0.2 Myr (lower panel). 


chemical model follows the detailed evolution of each SPH 
particle and the abundances of different elemental species 
are consistently derived at each timestep using the life¬ 
time function (Padovani & Matteucci 1993) and metallicity- 
dependant stellar yields from core-collapse SN, AGB stars, 
SNIa and pair instability SN. 

Galaxies are recognized as gravitationally-bound groups 
of at least 32 total (DM+gas+star) particles by running a 
friends-of-friends (FOF) algorithm. Of the galaxies identi¬ 
fied in the simulation, in our calculations we only use ‘well- 
resolved’ galaxies that contain at least 10 star particles, cor¬ 
responding to more than 145 total particles at any redshift. 
As a result, the number of galaxies is 1000 at z ~ 7 and 
1863 at z ~ 6. 

For each galaxy, we compute the intrinsic spectral 
energy distribution (SED) summing the contribution of all 
its stellar populations by means of the Starburst99 code 
(Leitherer et al. 1999; Vazquez & Leitherer 2005). In doing 
so, we consider the ages and metallicities of each stellar 
population as predicted by the numerical simulation. Then 
we convert the luminosity at A rea = 1500 A into absolute 
UV magnitudes, Muv- Given the UV magnitudes of the 
galaxies in the observational sample, in what follows we 
will focus only on the brightest galaxies in the simulation 
box, selecting the ones with Muv < —18 (102 objects at 
z ~ 7 and 225 at z ~ 6). The metallicities and stellar 
ages of these galaxies increase with the stellar mass, with 
5 x 10~ 3 ZQ < Z < 0.1 Zq and 40 Myr < t age < 200 Myr for 
7.5 < LogMstar/M® < 9.5. In addition, their star formation 
rates are < SOMg/yr, hence they can be considered as 


“normal” galaxies. 

The dust content of a galaxy depends on its star formation 
and metal enrichment histories. Since dust evolution is not 
computed within the simulation, our approach was to recon¬ 
struct the star formation history of each simulated galaxy 
at z = 6 and z — 7 along its hierarchical assembly and then 
use this as an input to the chemical evolution model with 
dust. This model has been first developed by Valiante et al. 
(2009), then applied to semi-analytical merger tree models 
of high-redshift QSOs by Valiante et al. (2011, 2012, 2014) 
and recently improved by de Bennassuti et al. (2014) to ex¬ 
plain the metal and dust content of the Milky Way and its 
progenitors. We refer the interested readers to these papers 
for more details. 

The equation describing the dust mass evolution is, 


Md(t) = —Z d (t)SFR(t) + Y d (t) 

Td T a cc 


( 2 ) 


where M d is the dust mass, Z d (t) = Md(t)/MisM(t) is the 
total dust abundance in the ISM, X c is the cold gas mass 
fraction and Td and r acc are the timescales for grain destruc¬ 
tion and grain growth and will be discussed below. 

The time-dependent term, Y d (t), is the dust production 
rate and it depends on the stellar IMF and on the adopted 
model for the stellar yields. We compute it as follows: 

rm up 

Y d (t) = / m d (m, Z) 4>(m) SFR(t — r m ) dm, (3) 

J m(t) 
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where md is the dust mass yield, which depends on the stel¬ 
lar mass and metallicity, $(m) is the stellar IMF, and the 
lower limit of integration, m(t), is the mass of a star with a 
lifetime r m = t. To be consistent with the assumptions made 
in the simulation, we adopt a Salpeter law in the mass range 
0.1 M@ < m < 100 MeQ. For stars with m < 8 Mg, we adopt 
metal yields from Van den Hoek & Groenewegen (1997) and 
dust yields by Zhukovska et al. (2008). For massive stars 
(12 Mg < m < 40 Mg) metal and dust yields have been 
taken from Woosley & Weaver (1995) and from Bianchi & 
Schneider (2007), including the effect of the reverse shock 
on dust survival in SN ejecta. For stars in the intermediate 
mass range 8Mg < m < 12 Mg, we interpolate between the 
AGB yields for the largest-mass progenitor and SN yields 
from the lowest-mass progenitor. Above 40 Mg, stars are 
assumed to collapse to black hole without contributing to 
the enrichment of the ISM. 

Finally, the last two terms in the right-hand side of 
eq. [2] represent the effects of dust destruction by interstellar 
shock waves and grain growth in the dense phase of the 
ISM. Here we assume that at any given time, a fixed fraction 
X c = 0.5 of ISM is in the cold phaseQ The timescale for grain 
destruction, ra, and grain growth, r a cc, are computed as in 
de Bennassuti et al. (2014). The latter timescale depends on 
the grain size distribution, and on the gas phase metallicity 
Z, temperature T mo i, and density n mo i, of the cold molecular 
phase where grain growth is most efficient, 



where we have assumed that grains which experience grain 
growth have a typical size of ~ 0.1pm (Hirashita et al. 
2014). For gas at solar metallicity, with n mo i = 10 3 cm -3 
and T mo i = 50 K, the accretion timescale is r aC c,o = 2Myr 
(Asano et al. 2013). Smaller grain sizes increase the grain 
surface area per unit dust mass, thus shortening the accre¬ 
tion timescale (Kuo & Hirashita 2012). Since the simulation 
has not the resolution to trace the temperature and density 
of the ISM, the values of n mo i and T mo i which best describe 
the galaxies under investigation are unconstrained. For this 
reason, in what follows we explore the impact of different 
values of r ac c,o on the resulting dust masses. 


4 RESULTS 

Figure |T] shows a comparison between the predicted dust 
masses of the simulated galaxies and the observations. We 
show the dust mass as a function of the stellar mass for all 
the simulated galaxies with intrinsic Muv < —18 at 2 = 6.33 
and s = 7.14. For each galaxy, the light grey point represents 
the mass of dust that is produced by stellar sources (AGB 
stars and SNe) after grain destruction and astration (eq. 0 


4 Since Pop III stars in the simulated galaxies are always largely 
subdominant, we neglect their contribution to chemical enrich¬ 
ment (Dayal et al. 2010). 

5 The value of X c , averaged over simulated galaxies with different 

masses, ranges between 0.38 and 0.56 (Vallini, Dayal & Ferrara 

2012 ). 


without grain growth), and the blue point is the total dust 
mass, including the effect of grain growth in the ISM. Here 
we have assumed a grain growth timescale r acc ,o = 2 Myr. 
Due to the lower metallicities, grain growth is not efficient 
at stellar masses M s t ar < 10 8 Mg. In these small galaxies, we 
predict that all the existing dust mass is entirely contributed 
by stellar sources. Galaxies with larger stellar masses have 
experienced a larger degree of chemical enrichment and the 
gas metallicities are high enough to activate efficient grain 
growth, although with a large scatter. Yet, due to the time 
required to enrich the ISM with metals, the bulk of the dust 
mass in the most massive galaxies is grown between z ~ 7 
and z ~ 6. In fact, while at z ~ 7 all the simulated galaxies 
have Mdust < 2 x 10 6 Mg, at z ~ 6 there are several galaxies 
with dust masses 10 7 M Q < Md us t < 10 8 Mg. Note that the 
observational data points, that we have reported in both 
panels for the sake of comparison, refer to galaxies which 
span a redshift range 6.5 < z < 7.5 (see Table 1). 

These results depend on the adopted stellar dust yield 
and grain growth timescale. In the lower panel of Figure[l]we 
show the effect of reducing the efficiency of grain destruction 
by the SN reverse shock: if all the grains survive the passage 
of the reverse shock, the effective SN dust yields is ~ 20 
times larger and the resulting dust masses contributed by 
stellar sources increase by a comparable factor. These SN 
yields bracket the observations of dust masses in SN and 
SN remnants obtained with the Herschel satellite (see Fig. 6 
in Schneider et al. 2014). The figure shows that although 
the total dust mass is larger, even with maximally efficient 
SN dust enrichment the predicted dust masses at z ~ 7 are 
Mdust < 10 7 Mg, too small to account for the observed dust 
mass in A1689-zDl. We conclude that in order to account 
for the existing dust mass in this galaxy, a shorter grain 
growth timescale is required. 

It is interesting to analyse the relative contribution of 
AGB stars, SN and grain growth as a function of redshift. 
Figure [2] shows the redshift evolution of the dust mass for 
galaxies with M s t ar > 10 9 Mg. We separate the contribu¬ 
tion from AGB stars, stellar sources, and the total mass 
of dust, including grain growth with an accretion timescale 
r ac c,o = 2Myr (upper panel) and 0.2 Myr (lower panel). As 
expected from Figure [T] grain growth provides the dominant 
contribution, exceeding the dust produced by stellar sources 
already at z < 10 — 12. Among the stellar sources, SN appear 
always dominant, but the average contribution of AGB stars 
can be as large as ~ 40%. This confirms that the contribu¬ 
tion of AGB stars to high redshift dust formation can not be 
neglected, especially for the galaxies currently targeted by 
observational searches (Valiante et al. 2009). Finally, in the 
lower panel we show that the dust mass detected in A1689- 
zDl requires very efficient grain growth, with a timescale 
T ac c ,o = 0.2 Myr, one order of magnitude shorter than re¬ 
quired to reproduce the observed dust-to-gas ratio in the 
Milky Way and in most local dwarf galaxies (de Bennassuti 
et al. 2014). This, in turn, implies that the cold atomic and 
molecular phases of the ISM, where grain growth is more effi¬ 
cient, must have an average density of ~ 10 4 5 cm -3 (see eq. 
Such a value is comparable to the molecular gas density in¬ 
ferred from CO excitation analyses of starburst galaxies at 
comparable (although slightly smaller) redshifts (see Carilli 
& Walter 2013). Morever, the molecular phase in these high- 
redshift star-forming galaxies may also be warmer, leading 
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to somewhat lower T aC c,o- Although plausible, our study sug¬ 
gests that these conditions must be exceptional, as if they 
were to apply to all galaxies a z > 6.5, current upper lim¬ 
its on the dust continuum emission for normal star forming 
galaxies at 6.5 < z < 7.5 would be exceeded. 


5 CONCLUSIONS 

In this work we have investigated dust enrichment of the ISM 
of “normal” galaxies, which form stars at rates between a 
few to a few tens of solar masses per year, at z > 6. We find 
that dusty galaxies are already formed at these high red- 
shifts, due to rapid enrichment by SN and AGB stars. How¬ 
ever, the dominant contribution to the dust mass in the most 
massive galaxies, with M star > 10 9 Mq, which are the galax¬ 
ies detected in the UV with present facilities, comes from 
grain growth. The efficiency of this process depends on the 
metallicity, density and temperature of the cold and dense 
phase of the ISM. While the upper limits on the dust contin¬ 
uum emission can be matched with an accretion timescale 
of T acc ,o ~ 2 Myr, comparable to the value adopted by de 
Bennassuti et al. (2014) to reproduce the dust-to-gas ratio 
of the Milky Way and local dwarfs, the observed dust mass 
of A1689-zDl at z = 7.5 can only be reproduced if more ex¬ 
treme conditions are adopted, with r ac c,o = 0.2 Myr, similar 
to the value assumed for QSO host galaxies (Valiante et al. 
2014). It is interesting to note that physical conditions in 
the ISM can lead to different values of r acc ,o even for galax¬ 
ies with comparable metallicity: two of the most metal-poor 
local dwarfs, SBS 0335-052 and IZwl8, have a metallicity 
of Z ~ 3% Zq but widely different dust masses (Hunt et 
al. 2014). The proposed interpretation is that the greater 
dust mass in SBS 0335-052 is due to the more efficient grain 
growth allowed by the high density ISM, observationally in¬ 
ferred to be almost 20 times higher than in IZwl8 (Schnei¬ 
der, Hunt & Valiante 2015). We conclude that, when the 
ISM is characterized by large enough densities, dust masses 
in high redshift galaxies may be more common than hith¬ 
erto believed and that the level of dust enrichment may not 
always be satisfactorily assessed from the gas metallicity. 
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